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ABSTRACT 

We present our catalog of U, B, V, and / stellar photometry of the central 64 deg 2 area of the 
Large Magellanic Cloud. Internal and external astrometric and photometric tests using existing 
optical photometry (U, B, and V from Massey's bright star catalog and I from the near-infrared sky 
survey DENIS) are used to confirm our observational uncertainty estimates. We fit stellar atmosphere 
models to the optical data to check the consistency of the photometry for individual stars across 
the passbands and to estimate the line-of-sight extinction. Finally, we use the estimated line-of-sight 
extinctions to produce an extinction map across the Large Magellanic Cloud, confirm the variation 
of extinction as a function of stellar population, and produce a simple geometrical model for the 
extinction as a function of stellar population. 

Subject headings: Magellanic Clouds — galaxies: photometry — galaxies: stellar content - 
dust, extinction — catalogs 



1. INTRODUCTION 

A galaxy's star formation history is encoded within its 
stellar populations. Outside of our own galaxy, the Mag- 
ellanic Clouds provide our most detailed view of galax- 
ies that are still vigorously forming stars. We present 
the catalog of stellar photometric data from the Magel- 
lanic Clouds Photometric Survey (MCPS) for the entire 
Large Magellanic Cloud (LMC) survey region (roughly 
8.5° x 7.5°, with the longer direction corresponding to 
the east- west axis). Our data are either deeper, cover a 
wider area, or include a larger number of filters (the in- 
clusion of U is particularly important for studies of dust 
and young stellar populations) than available optical cat- 
alogs. 

In addition to describing and providing the catalog, 
we construct and analyze extinction maps of t he LMC. 
As w e demonstrated fo r a portion of the LM C ijZaritskvl 
11999ft and for the SMC (|Zaritskv et alJl2002D . the extinc- 
tion properties in the Clouds vary both spatially and as a 
function of stellar population. Therefore, for many scien- 
tific purposes the catalog alone is insufficient, one must 
correct the observed magnitudes and colors for a com- 
plex extinction pattern. We describe the MCPS in §2, 
use the photometry to generate extinction maps of the 
LMC for two different stellar populations in §3, model 
the relative distributions of young stars, old stars, and 
dust in §4, and summarize in §5. 

2. THE DATA 

The data come from Magella nic Cloud Photometri c 
Survey, originally described by iZaritskv et al.l |l997). 
Using the Las Campan as Swope telescope ( lm) and the 
Great Circle Camera l|Zaritskv et al.lll996|) with a 2K 



CCD, we obtained drift-scan images for both Magellanic 
Clouds in Johnson U,B,V, and Gunn i. The effective 
exposure time is between 3.8 and 5.2 min for LMC scans 
and the pixel scale is 0.7 arcsec pixel -1 . Typical see- 
ing is ~1.5 arcsec and scans with seeing worse than ~ 
2.5 arcsec are not accepted. Magnitudes are p laced on 
the Johnso n-Kron-Cousins photometric system ijLandoltl 
1198^11991 . Scan images from observing runs from Oc- 
tober 1995 to December 1999 are included in this catalog. 
Additional pointed observations were performed in De- 
cember 2001 to fill in for unacceptable data from earlier 
observations. The data are reduced u sing a pipeline that 
utilizes DAOPHOT II f|Stetsonlll98l and IRAF 1 . Only 
stars with both B and V detections are included in the 
final catalog. 

Details of the reduction procedure, construction of 
the catalog, and qu ality assurance are presented by 
IZaritskv et all l)2002fl . Because the SMC and LMC data 
were taken concurrently, the quality of the two datasets 
is identical. We present only a cursory description of 
the procedure here. The pipeline for reducing individ- 
ual scans is a fairly standard application of DAOPHOT. 
Each of the 76 scans (x 4 for the four filters), which are 
11000 pixels long and 2011 pixels wide, is divided into 11 
by 2 subscans that are roughly 1100 by 1100 pixels, with 
overlap between the subscans that enables us to compare 
the results from the independent photometric reductions. 
The result of the reduction pipeline is a catalog of instru- 
mental photometry for each detected star in each filter 
and its right ascension and declination. We derive the 
astrometric solution from a comparison to stars in the 

1 IRAF is distributed by the National Optical Observatories, 
which are operated by AURA Inc., under contract to the NSF 




Fig. 2. — The stellar density map of the LMC for the young 
population (V < 18.5 and -0.2 < B — V < 0.3) in the left panel 
and an old population (V < 19.1 and 0.7 < B — V < 1.5) in the 
right panel. North is up and East is to the left. The image shows 
the entire ~ 8.5° X 7.5° survey region. 



Fig. 1. — A section of the red clump stellar density map for the 
LMC. The left panel contains the I-band photometry with no cor- 
rections, the right panel contains the photometry after correction. 



Magel lanic Catalogue of Stars (MACS; iTucholke et all 
(1996)), whose coordinates are on the FK5 system. 

We match the instrumental magnitudes of stars in dif- 
ferent filters using a positional match that associates the 
nearest star on the sky within an aperture that is 3 times 
either the positional rms of the astrometric standards in 
that subscan or 1.2 arcsec, whichever is larger. The V 
frame is used as the reference and only stars that have a 
match in the B frame are retained for the final catalog. 
In crowded areas it is possible that the "nearest" star 
in one filter is not the correct match to the V reference 
because of the uncertainties in the astrometric solution. 
We see some evidence of this problem when comparing 
to other data and when fitting atmospheric models (we 
find stars with highly anomalous colors). Except near 
the faint limit of the catalog or in extremely crowded re- 
gions, this issue appears to be a minor problem. These 
errors can be estimated reliably using artificial star sim- 
ulations. 

One significant difference bewteen the photometric so- 
lution of the LMC and SMC data, is that for the LMC 
data we do not automate the adjustment of the scan-to- 
scan photometry. We find that the photometry among 
scans is generally consistent, and only a few manual ad- 
justment of scan zeropoints were necessary. These were 
implemented to remove obvious discrepancies observed 
in a map of the density of red clump stars (for an exam- 
ple see Figure 0. The adjustment of specific individual 
scans avoids the potential for systematic drift present in 
an automated procedure. 

We present the catalog of astrometry and photometry 
for 24,107,004 stars as an ASCII table (see Table 1 for 
a sample). Columns 1 and 2 contain the right ascen- 
sion and declination (J2000.0) for each star. Columns 
3-10 contain the pairings of magnitudes and uncertain- 
ties for U, B, V, and I magnitudes. The last column 
is a quality flag that is described in Table 2, based on 
comparisons to Massey's (2001) catalog and our fitting 
of the spectral energy distribution (§3). A V band stel- 
lar density map of the LMC constructed from the catalog 
using stars with 16 < V < 20 is shown in Figure El The 
digital catalogs allow one to make analogous images for a 



variety of populations (for an example of the distribution 
of young and old s tars in the LMC, as done for the SMC 
bv lZaritskv et all l)2000|) . see Figure 2). 

The magnitude limit of the survey varies as a function 
of stellar crowding. We find little visible evidence for in- 
completeness at V < 20 (Figure EJ, but the scan edges 
and different scan sensitivities become visible when plot- 
ting the stellar surface density for stars with 20 < V < 21 
(Figure UJ. The U and / data are incomplete at brighter 
magnitudes than the B and V data. The U and / pho- 
tometry, even in sparse areas, is severely incomplete be- 
low U — 21.5 and I = 22 (comparable limits in the two 
other bands are B = 23.5 and V = 23). Any statistical 
analysis of this catalog fainter than V < 20 requires ar- 
tificial star tests to determine incompleteness, which is 
becoming significant at these magnitudes. 

2.1. Photometric and Astrometric Accuracy 

We presented extensive tests of the photomet ry and 
astrometry of the MCPS in lZaritskv et all l|2002j) . How- 
ever, we revisit some of those here for completeness. One 
particularly useful internal test is based on a map of the 
red clump mean luminosity. We calculate the mean mag- 
nitude of the red clump in 50 x 50 arcsec boxes over 
the survey region. Although large-scale variations in the 
mean magnitude may truly exists (for example due to a 
tilt of the LMC relative to a cons tant-distance surface; 
see Ivan der Marel fc Cioni I l)2001j) for a demonstration 
of an analogous effects with the red giant branch tip 
magnitude), any localized variation, in particular one 
that traces scan or subscan boundaries, reveals a prob- 
lem region. In Figure we show the maps of the de- 
viations from the mean red clump magnitude, for each 
filter. There are various important features in these pan- 
els. First, all four panels show increased noise toward the 
edges because there are fewer clump stars at large radii 
from the LMC and contamination for foreground Galac- 
tic stars is proportionally greater. Second, the U and / 
band frames show more structure that follows scan and 
subscan edges. Although we interactively corrected the 
most egregious of these (§2), we did not correct subscans 
within the bar region nor did we apply correction in cases 
where the correction needed was ambiguous. Recall that 
the U and I images are not as deep as the B and V 
images, so that some of the differences could be due to 
senstivity and completeness variations among scans (the 
clump stars are near the U magnitude limit). Third, 
there are large spatial scale variations that trace the bar, 
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TABLE 1 

The MCPS Catalog for the Large Magellanic Cloud a 



RA 


Dec 


U 


a v 


B 




V 








Flag 



4.490203 


—72.22987 


.000 


.000 


20.740 


.060 


19.746 


.054 


18.535 


.063 





4.490215 


-72.19868 


.000 


.000 


22.694 


.271 


22.788 


.310 


.000 


.000 





4.490216 


-72.28524 


.000 


.000 


21.646 


.099 


20.894 


.071 


.000 


.000 





4.490218 


-72.37505 


.000 


.000 


21.455 


.107 


21.050 


.086 


20.795 


.167 





4.490222 


-72.29113 


.000 


.000 


21.379 


.076 


21.298 


.106 


.000 


.000 





4.490228 


-72.09232 


.000 


.000 


21.668 


.104 


20.700 


.067 


20.114 


.087 





4.490231 


-72.09853 


21.724 


.410 


21.221 


.078 


20.740 


.062 


20.015 


.099 





4.490231 


-72.17493 


.000 


.000 


22.277 


.202 


22.220 


.213 


21.657 


.271 





4.490236 


-72.18102 


20.326 


.168 


20.013 


.070 


19.134 


.078 


18.055 


.073 


10 


4.490241 


-72.09316 


21.588 


.384 


22.491 


.231 


22.009 


.176 


.000 


.000 





4.490252 


-72.13191 


.000 


.000 


23.283 


.378 


22.997 


.356 


.000 


.000 





4.490256 


-72.13083 


20.287 


.140 


20.379 


.061 


19.989 


.055 


.000 


.000 





4.490261 


-72.20712 


.000 


.000 


22.687 


.318 


23.412 


.495 


.000 


.000 





4.490273 


-72.18983 


.000 


.000 


22.368 


.201 


21.895 


.135 


21.320 


.271 





4.490280 


-72.40957 


.000 


.000 


21.811 


.120 


21.296 


.113 


.000 


.000 






A The complete version of this table is in the electronic edition of the Journal. The printed 
edition contains only a sample. 



TABLE 2 
MCPS Catalog Flags a 



Description 



Value 



Replaced with Massey's (2001) photometry & astrometry 1 
Colors successfully fit with stellar atmosphere model 10 
Colors poorly fit with stellar atmosphere model 20 

A The flag values are additive. For example, a star with a flag value 
of 11 has had its photometry replaced with that of Massey's (2001) 
catalog and is fit well with the atmospheric models. 





and repeat in the various colors. These are real varia- 
tions in the mean magnitude of clump stars, presumably 
either due to physical effects (geometry, population dif- 
ferences, or variable extinction). We conclude that with 
the exception of some regions in the U and / band data, 
the photometry is self consistent among scans. The worst 
vari ations seen in the U b and correspond to ~ 0.1 mag. 

In Zarits kv et alJ l|2002D we compared our photometry 
to various existing catalogs. The quality of the photom- 
etry presented here is identical since the LMC and SMC 
data were taken concurrently. H owever there is one issue 
that requires further comment. iMasseyj (|20Q2f) has pro- 
duced a photometric catalog ([/, J5, V and R) of bright 
stars in the LMC. In particular, because of his interest 
in the upper main sequence, he pays acute attention to 
obtaining accurate t/-band photometry for young LMC 
stars, which is often difficult because of the lack of blue 
ill — B < 0.0) standards. He found that the combination 
of filter, CCD passbands, and the gravity sensitivity of 
the Balmer jump produced rather complex color trans- 
formation between his [/-band observations and the stan- 
dard Johnson system. Specifically, he required a second 
order color-term to correct his magnitudes. 

We had previously used his uncorrected photometry to 
correct our [/-band photometry (a slightly linear color 
term was added to place our photometry on his system, 
which was more extensively calibrated than our [/-band 




Fig. 3. — The stellar density map of the LMC for stars with 
16 < V < 20. North is up and East is to the left. The image 
shows the entire ~ 8.5° X 7.5° survey region. The nearly-horizontal 
white lines are caused by slight gaps between scans, which occurred 
if either the V or B scan was slightly misaligned. The irregular 
patchiness is due to extinction that is confirmed in §3. 



system; iZaritskv et al.l l|2002|) '). However, between that 
correction and his publication of the catalog, he identi- 
fied the additional complication in the calibration of the 
[/-band data mentioned above. Because the effect ap- 
pears to arise from the application of a calibration based 
on dwarf stars to the supergiants at these magnitudes 
in the LMC, our expectation is that our published pho- 
tometry is correct for dwarfs, but that the second order 
color term is required for supergiants. Because most of 
the stars in the catalog are not blue supergiants, our 
general results should be unaffected. This conclusion is 
supported indirectly by the results of our modeling of the 
SMC population, where we found go od agreement with 
the o verall photometry and models ({Harris fc Zaritskv! 
2003). Nevertheless, detailed photometry of the upper 
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Fig. 4. — Stellar density map of stars with 20 < V < 21. Various 
signs of incompleteness and double detections that were unmatched 
are evident. The scan-to-scan variations reflect varying complete- 
ness limits due to seeing variations. Bright regions of overlap be- 
tween frames show the incidence of false detections. The dark 
patches along the bar and the region of 30 Doradus show crowding 
compleness limits. 
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Fig. 5. — Comparison of mean red clump magnitudes in U 
(top left), B (top right), V (bottom left), and I (bottom right). 
Greyscale ranges from —0.3 to 0.3 mag. The entire survey region 
is shown. 



main sequen ce is more uncertain t han quoted in our orig- 
inal catalog ijZaritskv et alJl2002|) . 

To examine this issue furt her, we match our photom- 
etry of LMC stars to that of lMassevI (|20Q2fl . We match 
stars brighter than V = 15 in Massey's catalog to stars 
in the MCPS by finding the star in the MCPS within 
7.5 arcsec that has the closest V magnitude, as was done 
for the SMC data. The selection of the search radius 
is a compromise between assuring that the search finds 
the corresponding star while moderating the incidence 
of false matches. Because the bright star catalog has a 
relatively low density on the sky, we use a large search 
radius. As expected, we find that our uncorrected U 
magnitudes have a strong color term relative to those of 
Massey (see Figure EJ. We can correct our magnitudes 
by roughly fitting a second order B — V color term, which 
produces the second panel in Figure El The correction to 
our U photometry is AU = -0.5(B - V - 0.7) 2 + 0.23 if 
the (B -V) <1 and AU = 0.05 otherwise. The panels 
on the right hand side of the Figure demonstrate that 
the correction does indeed decrease the scatter between 
the two sets of measurements, but the systematic error is 
still as large as 0.1 mag at certain values oi B — V. Given 
the uncertainty of which stars the correction should be 
applied to due to the gravity sensitivity, we have cho- 
sen to not apply the correction to our data as presented 
in the catalog. However, we will present results in the 
discussion of derived extinctions with and without the 
correction. We find that whether the correction is ap- 
plied or not does not affect our results. 

Stars in our catalog that are brighter than 13.5 in B 
or V are prone to s ubstantial photometric uncertainty 
ijZaritskv et alJl2002|) . In the SMC catalog, we replaced 
the photometry and ast rometry f or sta rs brighter than 
this limit with those of iMassevI l|2002fl . However, the 
Massey catalog for the LMC covers only about a quar- 
ter of the area that ours covers. Therefore, most of our 




B-V A Color 



Fig. 6. — The comparison of magnitudes as a function of stellar 
color between the Massey survey and the MCPS. The upper panel 
shows the results from the original U band photometry, the second 
panel shows the corrected U band photometry. 



stars brighter than 13.5 in B or V have suspect photom- 
etry and astrometry. We have, when possible, replaced 
our photometry with Massey's, applying the mean pho- 
tometric offsets found between the MCPS and Massey's 
catalog, to Massey's data. To indicate which stars do 
have corrected photometry the catalog includes a flag 
(see Table 2 for description of the quality flags). The 
photometry and astrometry of 10912 stars are corrected 
using Massey's catalog. 

Finally, we compare our /-band photometry to that 
in the DENIS catalog. Although the DENIS catalog is 
prim arily an IR catalo g, it contains an / band channel 
and lCioni et alJ l)2000ft present a point-source catalog in 
the regions of the Magellanic Clouds. We use a search 
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Fig. 7. — The comparison of DENIS / band catalog with our 
data. The panel in the left illustrates the astrometric differences 
for matched stars with /-band magnitudes. The panel in the right 
shows the distribution of differences in I magnitudes between the 
two surveys. 



aperture of 3.5 arcsec for matches (smaller than that for 
comparison to the bright star catalog because of the in- 
creased stellar density). The distribution of astrometric 
and photometric differences for matched stars are plot- 
ted in Figure H In agreement with our previous results, 
we find that the astrometric accuracy is subpixel for the 
majority of the matches. The mean difference is 0.45 arc- 
sec. For matched stars whose magnitude agree to within 
1 magnitude, the mean photometric difference between 
the two surveys is 0.01 mag. 

3. EXTINCTION PROPERTIES 

A comparison between stellar atmosphere models and 
observed colors can b e used to infer the extinction to- 
ward individual stars l)Orebel fc I^ober t£lll995t IZa ritskv 
1999). We use the models of iLeieune et all i|1997fl and 
our U, B, V, and I photometry to measure the effec- 
tive temperature, T E , of the star and the lin e-of-sight 
extinc tion, Ay, in the manner described by IZaritskv! 
(1999). We a dopt a standard Galactic extinction curve 
(Senile! Il977j) . which is accepta ble for the LMC oy er 
the optical wavelength region fsee lOordon eTaTl IpTOl V 
The model fitting is least degenerate between T E and 
Ay for stars with derived temperatures in the ranges 
5500 K < T E < 6500 K and 12, 000 K < T E < 45000 K. 
Therefore, we construct Ay maps of the LMC from the 
line-of-sight Ay measurements to the set of "cool" stars 
(5500 K < T E < 6500 K) and the set of "hot" stars 
(12000 K <T E < 45000 K) with good quality photom- 
etry in all four filters (<j(j < 0.4, a E < 0.2, cry < 0.2, 
a i < 0.2) and good model fits (x„ < 3). In addition 
to these criteria, we imposed a reddening-independent 
magnitude cut {V < 19.0 + 3.2 x (B - V)). We caution 
however, that although the magnitude cut is reddening- 
independent, there is a bias in the catalog against highly 
extincted stars simply because the catalog itself is magni- 
tude limited. Stars which satisfy the photometric criteria 
and have an acceptable model fit (xt < 3) have a qual- 



ity flag of 10 and those that do not have an acceptable 
model have a quality flag of 20 in our catalog (Table ■ 

In Figure |H1 we show the spatial distribution of line- 
of-sight extinction derived from both the hot and cold 
stellar populations 2 . Because the recovery of Ay is quite 
sensitive to color, subtle differences in the scan photome- 
try are highlighted in the extinction maps. For example, 
a set of small photometric differences (0.03 mag in oppo- 
site senses in B and /, so that B — I has changed by 0.06 
mag) creates an extinction discontinuity of the magni- 
tude observed at many scan edges in the hot population 
map of Figure |H1 

The principal coherent extinction structure within the 
LMC is the increase in extinction in the hot star popula- 
tion along the north-east ridge of the LMC bar. This 
structure is also visible in the map from the colder 
stars. As expected, the increased extinction correlates 
with sites of star formation that can be seen in the 
stellar density images (Figure EJ), in lOO ^m images 
(iWheelock et al.lll99ll) . and in CO maps l|Fukui et alJ 
1999). 

The Ay histograms of the two populations are shown 
in Figure ED As we found for a small region of the LMC 
l)ZaritskvlTl999]) . the mean extinction is lower for the 
cooler populations (average Ay's of 0.43 mag vs 0.55 mag 
for the cold vs. hot population, respectively). The bi- 
modal distribution of extinction values among the cooler 
stars is characteristic of a geometry where the stars are 
distributed in front and behind a thinner mi d-plane dust 
plane (see Figure 12 of lHarris et~aT] l|1997^ . The lack 
of such bimodaility in the Ay distribution of the hotter 
stars suggests that those are distributed within the dust 
layer. This geometric model is further supported by the 
correspondence between the Ay values of the maximum 
in the hot star Ay distribution and the local minimum 
in the cold star Ay distribution. 

3.1. A Simple Model 

To explore the dependence of the Ay distribution on 
the relative distributions of the dust and stars, we de- 
velop a simple geometrical model. We construct a galaxy 
using three components: a "cold" star component that 
matches the projected spatial distribution of such stars, 
a "hot" star component that matches the projected spa- 
tial distribution of such stars, and a dust layer that corre- 
sponds to the dust map developed from the "hot" stars 
plus additional components for a diffuse, homogeneous 
internal dust layer and Galactic foreground, which is es- 
timated to be minimal {Ay ~ 0.05). The two dust com- 
ponents associated with the LMC, a clumped component 
seen in the extinction map and a possible smooth com- 
ponent that is not as easily detected, are defined to have 
the same vertical scale height. We then randomly place 
"cold" stars along the line of sight, but spatially con- 
strain them to follow the observed distribution by plac- 
ing them at the same (a, 5) as each observed "cold" star. 
Each star is reddended according to the local value of 
Ay (taken to reflect the midplane value of extinction as 
derived from the hot stars) and its position relative to 
the dust layer (for example, stars in front of the disk 
layer are extincted only by the Galactic foreground com- 

2 The extinction catalog is available for querv through the MCPS 
home page at http://ngala.as.arizona.edu/dennis/mcurvey.html 
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Fig. 8. — The spatial map of extinction values as derived for both the cooler, older stars (5500 K < Te < 6500 K ; right panel) and the 
hotter, younger stars (12000 K < Te < 45000 K; left panel) for the entire survey region. The small localized circular regions of apparent 
high extinction are globular clusters, which have some anomalous photometry because of their high stellar densities. Sharp spatial variations 
are due to scan-to-scan photometry differences of a few hundredths of a magnitude. 
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Fig. 9. — The distribution of extinction values for the two stellar 
populations and a simple model. We show the distribution of Ay 
for the hotter stellar population (upper panel) and for the colder 
stellar population (lower panel), where the solid and dotted line 
represent the unadjusted and adjusted (7-band photometry. The 
dashed line in the lower panel represents the expected distribution 
of Ay for an older stellar populations derived from the simple 
model detailed in the text. The sharp rise at zero extinction is 
due to observational uncertainties that scatter measurements to 
negative Ay, which we set to Ay = 0. The model shows the same 
rise at Ay = 0. 



ponent, while stars behind the dust layer are extinctcd 
by twice the measured value of Ay plus the hypothesized 
diffuse layer. 

We explore a range of relative thicknesses of the star 
and dust layers, and of the optical depth of the diffuse 
layer. A model that reproduces the salient features of 
the distribution of Ay for the cold population is shown 
in Figure [!|] The particular model shown has a scale- 



height for the cold stars that is ten times that of the 
dust, a Gaussian distributed error in the measured Ay 
of 0.125 mag, an optical depth for the diffuse layer, r/, 
of 0.28 (if the optical depth as estimated by the young 
stars in a region is < 0.28, we set the value to 0.28), 
and a foreground extinction optical depth of 0.05, which 
sets the position of the peak of the low- Ay population. 
The model excellently reproduces the Ay distribution 
for Ay < 0.7 and underestimates Ay for larger values. 
It may therefore be the case that there are pockets of 
higher optical depth in the midplane than modeled or a 
radial dependence in the extinction of the diffuse com- 
ponent that perhaps increases toward the LMC center, 
but adding those in an ad hoc manner to the models in- 
troduces far too many free parameters. Because of this 
systematic issue, we do not provide "best-fit" values for 
the parameters of our model, but instead simply show 
one satisfactory example. 

There is some degeneracy between the various param- 
eters, but large deviations from the values quoted above 
produce qualitatively inferior fits. For example, modi- 
fying Tf to correspond to Ay = 0.4 produces a highly 
inferior fit (Figure If OJI and decreasing the stellar scale- 
height by a factor of two fills in the valley between the 
two peaks of extinction values. An interesting variant of 
the model is to presume that Tf has a radial dependence. 
If Tf decreases exponentially with increasing radius, we 
find that for small values of the scalelength (for example, 
a 1 kpc scalelength for the model shown in Figure IT7)|) . 
the second peak is diluted, due to the spread in midplane 
extinction values in such a model, but that we do better 
at reproducing the high end tail of the extinction distri- 
bution because we can drive the extinction at the center 
of disk high and not affect large numbers of stars. Of 
course there are other solutions to the high-end extinc- 
tion distribution, so this agreement alone does not argue 
for the model with exponentially declining Tf. Models 
with large scalelengths approach the uniform Tf model, 
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Fig. 10. — A comparison of the distribution of extinction val- 
ues for the older stellar population and models. The histogram 
shows the data, the solid curve shows the model shown in Figure 
ITfl the short-dashed line represents the Tf = 0.4 model, the long- 
dashed line represents the model with a stellar scaleheight that is 
5 times that of the dust layer, the dot-dahsed line represents the 
exponentially declining Tt model. 



24 million stars in the 8.5° x 7.5° survey area centered on 
the Large Magellanic Cloud. The catalog contains posi- 
tions (right ascension and declination in J2000 coordi- 
nates) and U, B, V, and / magnitudes and uncertainties 
in the Johnson-Kron-Cousins photometric system mea- 
sured from our drift scan images. 

Using this catalog, we have constructed extinction 
maps for two stellar populations in the LMC. We find 1) 
that this dust is highly localized near the younger, hot- 
ter stars, and in particular toward regions immediately 
east and northeast of the center of the LMC, 2) that 
aside from these regions of higher extinction, there is no 
discernible global pattern, 3) that on average the extinc- 
tion toward the younger, hotter stars is only about 0.1 
mag larger, but that the distributions of extinctions are 
entirely different, 4) that the distribution of extinctions 
along lines-of-sight toward the older stars is bimodal, and 
5) that the bimodial Ay distribution is easily modeled 
as stars in front and behind a thinner dust layer. The 
two external galaxies for which we now have highly de- 
tailed maps of extinction as a function of stellar popu- 
lation both show significant differences in the extinction 
toward those populations. This difference, or at least the 
potential for this difference, should be considered when 
correcting the photometry of other galaxies for internal 
extinction. 



while models with smaller scalelengths further dilute the 
second peak. 

While the distribution of dust is undoubtedly more 
complicated than the description adopted in these mod- 
els, the models are able to reproduce the peak at Ay — 0, 
the peak at Ay ~ 0.05, the bimodal peak distribution, 
the position of the second peak, and the tail toward high 
Ay values with fairly minimal model assumptions. This 
suggests that as a global average the model is quali- 
tatively correct. Even so, it is a much more compli- 
cated description of the internal absorption in galaxies 
than generally ado pted. Some studies that focus on 
internal extinction (|Witt. Thronson. fc Capuanol 119921: 
iMisselt et al.ll200ll) explicitly deal with the distributions 
of stars and dust, but the general default correction is 
based on a foreground sheet assumption and an effective 
extinction curve. Because of the partial correlation of 
dust with star formation, this effective extinction curve 
is likely to be highly complex, dependent on geometry, 
and sensitive to the evolutionary state of the system. 
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4. SUMMARY 

We have conducted a broad-band photometric survey 
of the Magellanic Clouds. We present the data for over 
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